Rigidity sensing explained by active matter theory 
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The magnitude of traction forces exerted by living animal cells on their environment is a mono- 
tonically increasing and approximately sigmoidal function of the stiffness of the external medium. 
This observation is rationalized using active matter theory: adaptation to substrate rigidity results 
from an interplay between passive elasticity and active contractility. 



Living cells respond to mechanical as well as biochem- 
ical cues. Rigidity sensing designates the web of complex 
mechanisms whereby a cell will adapt, as a function of 
the elastic modulus of its environment, aspects of its phe- 
notype as diverse as its motility, gene expression, prolif- 
eration, or fate after differentiation [l|, 0| . 

The traction forces that a cell exerts on a flat, elas- 
tic plate depend in a non-trivial way on the extracellular 
stiffness [sj: roughly proportional to the elastic modulus 
in a softer environment, they saturate to a finite value for 
stiffer substrates [I]- Similar force- rigidity data are ob- 
tained whether forces are measured locally [1] or globally 
[^-0], for integrin-mediated 0, |3| as well as cadherin- 
mediated adhesion [1, |^, and even when the traction 
forces are exerted by assemblies of cells in a monolayer 
epithelium [l^. Single-cell rlieology assays @, show 
that the response of a cell to a sudden change of substrate 
rigidity is faster than the data acquisition rate allows to 
detect. These observations call for a simple, generic ex- 
planation, valid on short time scales where cell signaling 
cannot operate. 

In the context of adhesion-dependent mechanosensing, 
a simple 'two-spring model' was introduced in pre- 
dicting that stiffer environments lead to stronger trac- 
tion forces. A 'three-spring model' was later proposed 
in order to explain the stiffness-dependent orientation 
of stress fibers in adherent cells [l^, where contractility 
modulates cytoskeletal stiffness via a phenomenological 
polarizability coefficient. In this Letter, we formulate and 
solve a simpler model derived from active matter theory, 
a generic description of living matter where the mechano- 
chemical transduction due to molecular motors (activity) 
plays a central role [13^, J^] . We obtain a (static) force- 
rigidity relationship that agrees well with experimental 
data. We give an expression of the (dynamic) loading 
rate, expected to be valid on time scales too short for cy- 
toskeletal remodeling and protein recruitment to occur 
(t-Cl02 s 



MODEL 

Active matter theory formulates constitutive equations 
that take into account viscoelasticity, activity, the po- 
lar nature of cytoskeletal biopolymers [l^, and respect 
the principles of linear irreversible thermodynamics. In 
one spatial dimension, we write the stress of an homoge- 
neous, elastic and contractile material as the sum of an 
elastic and an active contribution: cr — cr°' -I- a a- The 
active stress > is proportional to A/i, the difference 
of chemical potential between products and reagents of 
the chemical reaction responsible for mechanochemical 
transduction (ATP hydrolysis): cta = C A^ t, where C is a 
material parameter of the cytoskeleton . For simplic- 
ity, we restrict our description of the cytoskeleton to the 
linear regime, and include: 

(i) elasticity, described by a linear spring of length 
lc{t) at time t, rest length Iq and spring constant 
kc; 

(ii) activity, modeled by an active force Fa = <Ja S 
across a section of area S. 

In the nonlinear regime, elastic moduli may also depend 
on activity, as hypothesized in • Assuming, as in [13] , 
that the rest length of elastic cytoskeletal structures is 
time-dependent due to motor activity is also beyond the 
scope of the linear regime that we consider here. 

The extra-cellular environment is represented by a lin- 
ear spring of length lext{t) at time t, rest length l^^^, and 
spring constant fccxt (see Fig. [T|). The traction force F{t) 
exerted by the cell on its environment reads: 

F(i)=A:e,t(/cxtW-0- (1) 

The sign is chosen as F > for contraction. Under usual 
experimental conditions, the total length ^tot of the sys- 
tem (cell -I- substrate) is constant: ^tot — lc{t) + ^cxt(0- 
The force balance equation reads: 

-Fit) + kc{lcit)^fc)+FA = 0. (2) 
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FIG. 1. Schematic representation of the model 

STATICS 

The amplitude of the traction force at equilibrium 
reads: 



Sat 



fccxt + kc 

When fcoxt 3> kc, the traction force saturates to i^sat 



(3) 



Fsat =FA + kc {It 



Ct-l'c), (4) 

the sum of Fa and of residual stresses kcihot — Icxt ~ ^c)- 
We expect the ensemble average of residual stresses to 
cancel: < Fsat >=< Fa >■ As long as fcoxt < ^c, 
the traction force is a linear function of fccxt: F^'^ ~ 

(^Sat/fcc) fcext- 

Experimentally, a wide range of rigidities fcext can be 
obtained when the substrate is a dense array of cylin- 
drical elastomer micropillars whose stiffness depends on 
their radius and height. Depending on the coating pro- 
tein, traction forces are transmitted through integrin- 
mediated adhesions (with fibronectin) [3, Q or through 
cadherin- mediated adhesions (with N-cadherin) Q. 
Fig. [2] show that, in both cases, experimental data are 
well fitted by the force-rigidity relation Eq. ([3]). Note 
that the fitted values of saturation forces i^sat were not 
observed: they correspond to values of /ccxt so large that 
the deflections of pillars would fall below the experimen- 
tal spatial resolution. 

The data encompasses three cell types, Madin-Darby 
Canine Kidney (MDCK) cells, 3T3 fibroblasts, and C2 
mouse myogenic cells. Cytoskeletal organisation and ad- 
hesive properties of the cells vary substantially according 
to the type of adhesions and the range of substrate rigid- 
ity. More diffuse cortical actin dominates when the envi- 
ronment is softer, wheras actomyosin bundles are prefer- 
entially formed at higher rigidities In all cases, our 
simple model captures the essence of the force-rigidity 
dependence, and sums up biological variation into two 
quantitative parameters, the asymptotic traction force 
Fgat and the cytoskeletal stiffness kc- The order of mag- 
nitude of the saturation force Fgat ~ 10 nN corresponds 
to an active stress of the order of 10^ Pa. We obtain 
this value upon neglecting possible residual stresses, and 
using OA Fsiit/S, where the section 5 of a micropillar 
is of order S ^ 1 /im^. The cytosketal rigidity is of the 
order of 10^ nN//im, corresponding to elastic moduli of 
the order of 10^ Pa, a value intermediate between moduli 
typical of cortical actin [Ts] and of stress fibers [llj (we 
used E ^ kc d/ S with d ~ 1 /im). 

When cells are treated with blebbistatin, an inhibitor 
of contractility, the value of Fsat is halved (Fig. [2]B). In- 
spection of Eq. (jH) suggests that traction forces remain 
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FIG. 2. We fit experimental data (average traction force per 
pillar, markers) with Eq. ((3| (lines). A: integrin-mediated 
adhesions (Figs. 2a and 2c of |3]). MDCK cells (black cir- 
cles): Fsat = 34 ± 4 nN, kc = 146 ± 30 nN//im (black 
solid line); fibroblasts (blue squares): Fsat = 19 ± 3 nN, 
kc = 119 ± 35 nN//im (blue dashed line). B: cadherin- 
mediated adhesions (Fig. 4a of [3|). Standard conditions 
(black circles): _Fsat = 29 ± 10 nN, kc = 209 ±112 nN/^im 
(black solid line); cells treated with blebbistatin (blue circles): 
Fsat = 15 ± 8 nN, kc = 183 ± 164 nN/^im (blue dashed line). 
Error bars correspond to 95% confidence levels. 



non-zero due to residual activity of myosin motors, as 
proposed in Q, or to non-zero residual stresses, or to 
a combination of both effects. We note that the value 
of kc is almost unchanged. However, assuming as in 
[T^ that activity modulates cytoskeletal stiffness via a 
polarizability coefficient a leads to a ratio F^'^/i^sat = 
^ext/ (fcext -I- (1 -|- a) kc) that depends through a on cy- 
toskeletal contractility. Experimental data [9j suggests 
that F°'^/Fsa.t is independent of the level of contractility, 
in agreement with our prediction, Eq. ([3]). 



DYNAMICS 

The loading rate exerted by T cells immediately af- 
ter receptor engagement with a model antigen-presenting 
cell was measured thanks to a biomembrane force set-up 
where micropipette aspiration controls the external rigid- 
ity, and found to be linear in fcoxt [20] . Motivated by this 
result, we turn to the dynamics of traction forces, and 
modify Eq. ^ by taking into account internal protein 
friction in a linear force- velocity relationship: 



FA{t)^Fs + ^ 



(5) 



where Fs is the stall force, and ^ is a friction coefficient 
[2l|. Eliminating other variables in Eq. ([2]), we obtain a 
differential equation for the traction force: 



dF 
"dF 



F 

T 



T 



(6) 



with a viscoelastic time t — ^/ (fccxt + kc)- Integration 
from an initial time tg with initial force F{t = Iq) = f ° 
gives F{t) = F<=i-)-(i^" - F^i) e-(*-*«)/^. For zero initial 
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force, we find that the initial loading rate is proportional 
to the substrate rigidity: 

-7- ^ fcoxt, (7) 

dt |(t-to)<ST C 

in agreement with [20 |. where the initial time to is set 
when pulling starts so that = 0. We checked that 
Eq. ^ still holds if we replace the linear force- velocity 
equation ^ by Hill's equation [22||. Since Fsat is a func- 
tion of activity, we predict that the loading-rate-rigidity 
data will be modified upon treatment with contractility 
agonists and antagonists. 

CONCLUSION AND OUTLOOK 

At low external rigidity, cell traction forces increase lin- 
early with the stiffness of the substrate. Their constant 
ratio was first interpreted as a displacement regulated by 
the cell We show that regulation is not necessary to 
explain the force-rigidity relationship. Within the frame- 
work of linear irreversible thermodynamics, we propose 
a minimal model whose consequences are consistent with 
available experimental data. We predict that both the 
saturation force Fsat, exerted for large stiffness, and the 
constant displacement i^sat/^d observed at low stiffness, 
depend upon the contractility level. Our description is 
relevant for several types of adhesive structures and cy- 
toskeletal organisation. Unlike [ll|, we ignore the dy- 
namics of adhesive contacts through which force is trans- 
mitted to the substrate. Other monotonically increasing 
functions of stiffness that depend on two parameters also 
fit the same experimental data: it is our hope that this 
work will foster further quantitative experiments to con- 
firm - or disprove - our model. 



To treat the dynamics of traction forces, we include in- 
ternal friction, and obtain an initial loading rate propor- 
tional to external stiffness, as observed experimentally 
[20| . In accord with single-cell rheology assays [1, 'rf, 
the loading rate d_F/dt responds instantaneously to vari- 
ations of fcoxt. 

The merit of our analysis is to show that the simplest 
equations of active matter dictated by symmetry and 
conservation laws are sufficient to describe a behavior 
which might seem to require a more elaborate regula- 
tion at first sight. This suggests that other features such 
as stress fiber diameter and equilibrium with the rest of 
the actin-myosin system could be described within the 
general framework of active gels. We hope that exten- 
sions of our model will allow to understand quantitatively 
how more complex cell processes depend on extracellu- 
lar rigidity. Including membrane elasticity and cortical 
tension in an appropriate geometry may explain why the 
initial loading rate exerted by T cells upon receptor en- 
gagement saturates for stiffer environements [20|. The 
dynamics of wetting of the microplate by the cell must 
be taken into account to describe traction forces exerted 
during cell spreading d, Q . Finally, biochemical signal- 
ing, protein recruitment and remodeling of adhesive and 
cytoskeletal structures act over longer time scales [H, [l^ 
and may enhance the mechanical effects described here. 
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